x = -1.5:0.01:1.5;
y = -1.5:0.01:1.5;
[X, Y] = meshgrid(x, y);
contour(X, Y, func(X, Y), 100);
v = logspace(-3, 0.7, 15);
hold on;
contour(X, Y, func(X, Y), v);
axis equal;
plot(1, 1, 'r*', 'MarkerSize', 20, 'LineWidth', 3);